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ABSTRACT 

We have examined consequences of strong tidal encounters between a neutron star and 
a normal star using SPH as a possible formation mechanism of isolated recycled pulsars in 
globular clusters. We have made a number of SPH simulations for close encounters between a 
main-sequence star of mass ranging from 0.2 to 0.7 M Q represented by an n=3/2 polytrope and 
a neutron star represented by a point mass. The outcomes of the first encounters are found to be 
dependent only on the dimensionless parameter rf = (m/(m + M)) 1 / 2 (r m i n /RMs) 3 ^ 2 ( m /MY 1 ^ 6 \ 
where m and M are the mass of the main-sequence star and the neutron star, respectively r m i n 
the minimum separation between two stars, and Rms the size of the main-sequence star. The 
material from the (at least partially) disrupted star forms a disk around the neutron star. If all 
material in the disk is to be acctreted onto the neutron star's surface, the mass of the disk is 
enough to spin up the neutron star to spin period of 1 ms. 

Subject headings: celestial mechanics - stellar dynamics - globular clusters: general - pulsars: 
general 
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1. INTRODUCTION 

The cores of globular clusters have high stellar densities. Recent studies showed that the physical 
interactions between stars are main driving force in determining the dynamical evolution after the core 
collapse (e.g., Goodman 1988). The physical interactions include direct collisions between normal stars, 
tidal captures between a normal star and a compact star (neutron star or white dwarf), and formation of 
binaries via three-body processes. 

Such interactions can lead to the formation of objects that are not common in low stellar environments. 
Indeed globular clusters reveal high abundance of X-ray binaries and short period pulsars. These systems 
are thought to be related to the physical interactions between stars. X-ray binaries can be formed either by 
tidal capture between a neutron star and a main-sequence star, or by the evolution of primordial binary. 

Another observational evidence for the high abundance of the physical interactions in clusters is the 
enhancement of short-period pulsars. A large number of short period pulsars have been detected in globular 
clusters. Unlike many pulsars in the solar neighborhood, cluster pulsars are known to be predominantly 
short period ones. Neutron stars are believed to be born with short spin periods (< 0.1 sec) and strong 
magnetic fields (10 12 — 10 13 G). As a neutron star evolves, its rotational energy is released through magnetic 
dipole radiation which causes deceleration of the spin. Since the spin-down time scales for pulsars with 
strong magnetic fields are of order 10 8 years while the production of neutron stars in globular cluster must 
have ceased long time ago, present day population of short period pulsars in clusters must be 'recycled' 
pulsars, which has been spun up by accretion of mass. 

Aside from the fact that the abundance of recycled pulsars in globular clusters relative to ordinary 
stars is far greater than that for solar neighborhood, the cluster pulsars are preferably isolated (as opposed 
to in a binary system). According to compilation of 558 pulsars by Taylor, Manchester, & Lyne (1993), 19 
out of 22 in clusters are isolated pulsars while the ratio is only 2 out of 11 in field pulsars. Therefore any 
theory for the recycled pulsars in globular clusters should explain both the high abundance and high ratio 
of isolated to binary pulsars. 

Proposed mechanisms for producing isolated recycled pulsars in dense stellar regions include close 
encounters between a neutron star and a main-sequence star leading to a complete or partial disruption 
of the main-sequence star (Krolik 1984, Lee 1992), and the three-body interactions which can liberate a 
neutron star in the binary through various ways (see Rappaport, Putney, & Verbunt 1989 for summary). In 
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the present paper, as a start of the study on the first mechanism, we examine the consequences of close 
encounters between a neutron start and a main-sequence star. 

During a close encounter involving a neutron star, the orbital energy of an incident normal star can 
be dissipated and stored in its interior by tidal force exerted by a neutron star. Depending on the amount 
of energy deposited to the stellar envelope, one of the following outcomes are possible: tidal capture of a 
normal star, or total or partial disruption of the normal star. Based on SPH simulations, Davies, Benz, 
& Hills (1992) find that encounters involving a main-sequence star of 0.8 M Q can leave a single object for 
fmin i5 1-75 Rms an d detached binaries for 1.7 Rms i5 r min i5 3.5 Rms, while encounters involving a red 
giant star leave a compact neutron star- white dwarf binary for r min < 1.8 Rue an d a detatched binary for 
1.8 Rrg ^ T'min iS 2.5 Rrg, where, r m ; n is the separation at periastron passage. In the present study, while 
Davies et al. concentrated on only one main-sequence star mass, we perform SPH simulations of encounters 
between a 1.4 M neutron star and a main sequence star with a mass in the range between 0.2 and 0.7 
M Q for various r m j n . This extension will bring a full comparison between simulation and theory, and will 
eventually provide more realistic information to the study of dynamical evolution of globular clusters. 

This paper is organized as follows. In §2, we describe our simulations of stellar encounters, and we 
discuss the consequences of the encounters based on the simulations in §3. Final section summarizes our 
major findings. 



2. SIMULATION 

We adopt an SPH code developed by Monaghan & Lee (1994) for the encounters between a massive 
black hole and a normal star which is modelled as an n=3/2 polytropic sphere. We use 9185 SPH particles 
to model a main-sequence star and treats the neutron star as a 1.4 M Q point mass that interacts with SPH 
particles only via gravitational forces. We have chosen main-sequence star mass m to be 0.2, 0.3, 0.5, and 
0.7 M© in our simulations, and the main-sequence stars are assumed to have a linear mass-radius relation 
(i.e. Rms/R-g = m/M Q ). For Rms of order R Q and relative velocity at infinity, Wqo, of order 10 kms -1 , 
the relative velocity v m i n at r m - m becomes 

2 G(M + m) 2 

"min = " + v oo (1) 

' min 

s ^M + m)^ (2) 
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where M is the mass of the neutron star. This means that Uqo is not an important parameter in determining 
the orbit of the star around the neutron star, and in our simulations, is alwasys the representative value 
in globular clusters, 10 kms -1 . Therefore we assume that the relative orbit is parabolic for simplicity. We 
made total of 16 simulations whose input parameters (m and r m in) are listed in Table 1. 

For the purpose of illustration of the outcome of tidal encounters, we have shown four snapshots of 
simulation m7p08 in Figure 1. Small dots are the SPH particles, and the curved line is the trajectory that a 
main-sequence star would follow if it were a point mass. In this figure, the coordinate center is the location 
of the neutron star. At t = 1.8 h (relative to the start of the simulation), the shape of the main-sequence 
star is already distorted like a Roche lobe which is pointing off the neutron star slightly. At t = 2.2 h, 
the SPH particles are spread out in two opposite directions so that the center of mass still reside near the 
parabolic orbit, and some of the particles from the main-sequence star start to spiral into the neutron 
star. At t = 3.3 h, a disk is formed around the neutron star and is still connected to the main body of 
the main-sequence star which has a long tail of particles with relatively smaller velocities. At t = 4.2 h, 
the distortion of main-sequence star has ceased and the main-sequence star is now composed of a sphere 
with high density and two spiral arms. The disk around the neutron star is expanding in the orbital plane 
because the particles in the disk has large eccentricities which make them spend much longer time at farther 
distances than the beginning of the accretion (t ~ 3.3 h) when most particles are near their perigees. 



3. ANALYSIS 



The main purpose of the simulations in this study is to observe how much mass of the main-sequence 
star can be transfered to the vicinity of (or bound to) the neutron star in the process of close encounter, 
as a function of m and r m j n (note that this mass may differ from the actual mass that accretes onto the 
neutron star surface and consequently spins up the star). Also, we will attempt to predict the future of the 
main-sequence star remnant, which can be parameterized by its size, r m i n at the second encounter, and the 
period of its orbit. The disk formed around the neutron star will be discussed in connection to the spin-up 
of the neutron star. 



3.1. Mass Fractions 
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To study separately the futures of the disk around the neutron star and of the main- 
sequence star remnant, it is necessary to judge whether each particle will be bound to 
the neutron star, remain as a part of the main-sequence star, or escape from the system. 
Here we employ the following criteria to distinguish among the above three possibilities: 

(A) if (Ens < 0) and (Ens — Ems < 0), the particle is bound to the neutron star; 

(B) elseif (Ems < 0), bound to the main-sequence star; 

(C) else, lost from the system. 

The relative energies Ens and Ems are defined as 

1 2 GM 

E NS = 2V 2 -— (3) 
Ems = i(v-(v>) 2 -^, (4) 

where v and r are the velocity and the distance of a particle relative to the neutron star and the brackets 
denote the averages taken over all SPH particles. Figure 2 is a reproduction of the last snapshot in Figure 
1 with three different symbols. Small dots, large dots, and x's are the particles that satisfy the conditions 
A, B, and C, respectively. Note that there are two clear separations between small dots and large dots and 
between large dots and x 's. The amounts of mass belonging to each criterion as a function of time for 
simulation m7p08 are shown in Figure 3. These masses are well defined because each line approaches to an 
asymptotic value rather quickly. 

However, the mass fractions by the above criteria for simulations m2p02 and m7p05 do not approach 
to asymptotic values, because their r m i n 's are so small and the main-sequence star experiences so strong 
tidal force that the main-sequence star sprials into the neutron star (or, more realistically, the latter spirals 
into the former) in a relatively short time instead of spending some time on the quite eccentric orbit as an 
independent object. Thus for these two simulations, it is impossible to distinguish the first encounter and 
the subsequent one. We do not include these two simulations in the following mass fraction anaylsis which 
is limited to the mass fractions of the first encounters only. 

The amounts of mass that satisfies criteria A, B, and C are plotted as a function of r m i„/i?MS m Figure 
4 for all simulations except simulations m2p02 and m7p06. Each symbol denotes different main-sequence 
star mass (see the caption). It is clear that at the same r m ; n /i?Afs, heavier main-sequence star loses less 
fraction of its mass. From the similarity of the curves for each membership, one may expect a parameter 
that soley determines the mass fractions. For a theoretical approach, we have first tried r\ of Press & 
Teukolsky (1977). The quantity r] measures the duration of periastron passage, relative to the hydrodynamic 
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timc of the star, and is defined as 

1/2/, . \ 3 / 2 
l / ' min \ 



m + Mj \Rms) 



(5) 



However, it is found that 77 does not describe the mass fractions in a unified way. Rather, we find that a 
small alternation to 77 can make the mass fraction be a function of only one parameter: As shown in Figure 
5, the mass fractions of all simulations allign on a single curve if the abscissa of Figure 4 is replaced with a 
new quantity 

^WTmJ teJ Km) ■ (6) 

Theoretical base of the above alternation is as follows: According to Press & Teukolsky (1977), when the 
relative orbit is assumed to be a parabola, the deposition of orbital energy to stellar oscillations of the star 
with mass Mi and the radius R\ due to the perturbing star with M 2 and R2 can be expressed by 

where I is the spherical harmonic index. The dimensionless function 7] (77) depends on the structure of 
the star, and the values of Ti(rj) for the star with a polytropic index n = 2/2, has been calculated by Lee 
& Ostriker (1986). Then one could expect that the amount of the mass detatched from a star during an 
encounter, Am, should be related to this AE. Furthermore, since the detatched mass is a function of 77' , 
AE is expected to be a function of 7/ somehow as well. Indeed we find that AE/m can be well expressed as 
a function of 77'. In Figure 6, AE/m's for three different main-sequence star masses are plotted over 77 and 
77'. This comparison clearly shows that AE/m can be expressed more uniformly in terms of 77' than just 77. 

It would be worth while to relate Am/m with AE/m. Such relationship can be found in Figure 7. 
Note that Am/m increases with AE/m linearly until Am/m ~ 0.1. After that, it increases quatratically, 
but the overall tendancy can be approximated to be still linear. 

The theoretical quantity AE can be obtained in our simulations by finding the difference of the orbital 
total energies (kinetic + potential) at the start and at the end of the simulation. However, since this 
calculation is meaningful only when there is almost no mass loss from the main-sequence star, we compare 
A_E's with large 77' only. In Figure 8, one can easily see that there is a good agreement between our 
simulations denoted by several kinds of symbols and the theory denoted by the curves. 



3.2. The Main-Sequence Star Remnant 
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All main-sequence stars in our simulations are captured by the neutron stars and thus will experience 
subsequent encounters. However, no numerical method has been known to be able to handle such long 
period of hydrodynamic simulations. For this reason, we have performed the simulations of first encounters 
only and try to predict the perigee distance at second encounter, r m ; nj2 , and the period between first and 
second perigees, P, from the last dumps of some of our simulations. 

To estimate r m i nj 2 and P, we take r and v averages for the particles that satisfy the criterion B in 
the previous subsection. Since the quantities r m i n .2 and P approach to asymptotic values less quickly than 
the mass fractions, we have extended simulations m7p08, m7pll, and m7pl5 until they show asymptotic 
values, which are listed in Table 2. 

Approximately, simulation m7pl5 has twice greater r m ; n than that of simulation m7p08, but its P is a 
thousand times longer implying that P has a strong dependency on r m j n . 

In the above three simulations, r m i nj2 is about 50 % bigger than r m i n . From this fact, the amount of 
mass stripped during the second encounter, Am2, might be exptected to be not as much as that for the 
first encounter, Am, but one should remind that the parameter that determines Am, rf', is a function 
of Rms too. Since the main-sequence stars are expanded during encounters (see Fig. 9), Rms will 
have been enlarged at the time of second encounter resulting smaller rj' and greater Am than would be 
with a constant Rms- Furthermore, as in case of simulation m7pl5, even considerable changes in the 
structure accompany. Also, the rotational angular momentum of the main-sequence star transfered from 
orbital angular momentum, which is another result of the tidal interaction, may play an important role in 
determining the amount of mass loss during subsequent encounters. Thus we conclude that the prediction 
for the fate of the main-sequence star remnant during subsequent encounters is very difficult to obtain from 
our simulations and that more polished numerical and/or method that can easily handle such a long period 
of hydrodynamic simulations is neccessary for that. 



3.3. The Disk 



Unlike the binary system where mass transfer to a disk from outside the Roche lobe of a companion is 
relatively stable, the formation of an accretion disk following a tidal disruption is so abrupt that too much 
material is put in the disk in a short amount of time. Consequently, the mass accretion rate from the disk 
onto the neutron star is much higher than the Eddington limit. Verbunt et al. (1987) estimated the amount 
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of accreted mass to be only (MEdd/ 'M)Md which is insufficient to lead to millisecond periods, where Msdd 
is the mass accretion rate corresponding to the Eddington limit luminosity, M the mass accretion rate onto 
the neutron star, and M d the disk mass. 

The spin period P s of the neutron star after the accretion is determined by total amount of accreted 
mass AMacc, and the the relationship between the final P s and AM acc can be written as 

/ p \ -4/3 

AM acc ~ 0.20/ I J M (8) 

( Phinney & Kulkarni 1994), where f > 1 is a factor that depends on the details of the accretion, and 
the moment of inertia of the neutron star is assumed to be 10 45 gcm 2 . Until the neuron star reaches its 
equilibrium period, / remains to be around 3 (Ghosh & Lamb 1992). Further accretion after the equilibrium 
period does not spin up the star since / has become much higher than 1. 

If the steady, thin disk analysis by Ghosh & Lamb (1992) with / ~ 3 is still applicable and if the 
estimation of mass transferred to the neutron star by Verbunt et al. (1987) is valid, the accretion followed 
by tidal disruption with a rate of 1000 MEdd may be able to accelerate P s to > 100 ms. 

On the other hand, if there is a way for all disk material to overcome the Eddington limit and to 
accrete onto the neutron star's surface, the accretion would be able to accelerate P s to ~ 1 ms. 

The reality will reside somewhere between the above two extreme cases, and the authors will pursue a 
study on this heavy accretion disk around a neuron star as a next step of the present study. 



4. DISCUSSION 

We have investigated the consequences of close encounters between a neutron star and a main-sequence 
star using the SPH method. We have found that the mass fraction stripped from the main-sequence star 
during the encounter, Am, is determined only by a single parameter i]' in the ranges of m and Voo used 
in our simulations. The orbital energy deposited in the stellar envelope per unit mass, AE/m seems 
responsible for the mass stripping because it also can be well approximated to be a function of rf. 

From simulations with relatively large rf , the core contraction of the main-sequence star remnant has 
been witnessed as well as the envelope expansion. The radiation from the contracted core will help the 
envelope expansion and the envelop may not be able to restore itself to the original place. 
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The mass of the disk formed around a neutron star after an encounter is as much as of order 0.1 M©. 
If all of this mass could be accreted onto the neutron star's surface, the accretion could accelerate the 
spin of the neutorn star down to 1 ms. However, this kind of heavy accretion disks are subject to high 
radiation pressure and considerable fraction of the disk is believed to be blown away, which is against the 
spin-up of the neutron star. The fate of these heavy accretion disks are uncertain and remains to be studied 
intensively, to explain the high abundance of isolated, recycled pulsars in globular clusters. 

We are grateful to Dr. M-G Park for useful conversation on the accretion disk theory. This research 
was supported in part by the Korea Research Foundation through non-directed research funds in 1993. 
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Table 1. Simulation Parameters 



Simulation 


M M s 
(M ) 


^*min 

(10 10 cm) 


V 


1 

V 


Simulation 


Mms 
(M ) 


^*min 

(10 10 cm) 


V 


V 


m2p02 


2 


2 


0.60 


0.44 


m5p07 


5 


7 


1.45 


1.22 


m2p03 


2 


3 


1.11 


0.80 


m5p09 


5 


9 


2.16 


1.78 


m2p04 


2 


4 


1.71 


1.23 


m5pll 


5 


11 


2.86 


2.41 


m2p05 


2 


5 


2.39 


1.72 


m7p05 


7 


5 


0.60 


0.53 


m2p06 


2 


6 


3.14 


2.27 


m7p08 


7 


8 


1.20 


1.07 


m3p04 


3 


4 


1.10 


0.85 


m7pll 


7 


11 


1.94 


1.73 


m5p06 


5 


6 


1.15 


0.97 


m7pl5 


7 


15 


3.09 


2.75 



Table 2. Orbital Elements of Selected Simulations 



Simulation r min r min , 2 a P 

(10 10 cm) (10 10 cm) (h) 



m7p08 8 10.9 23 

m7pll 11 18.1 210 

m7pl5 15 22.4 28000 



a Pcrigcc distance at the second encounter 
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Figure not included because of the large file size 

You may get this figure at 

http: // mercury, es.pusan.ac.kr/papers. html 



Fig. 1. — Four snapshots of simulation m7p08 with Mms — 0.7 M© and r mm = 7 x 10 10 cm. The hyperbola 
is the trajectory that the main-sequence star would follow if it were a point mass. All snap shots are in 
ncutron-star-relative coordinates such that the neutron star represented by a cross resides always at (0,0). 
The axes are in units of Rms- 
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Figure not included because of the large file size 
You may get this figure at 
http://mercury.es.pusan.ac.kr/papers.html 



Fig. 2. — Reproduction of the last snapshot in Figure 1. Small dots are the particles bound to the neutron 
star, large dots those bound to the main-sequence star, and x those unbound to the system. For membership 
criteria, see text. The axes are in units of Rms- 
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Fig. 3. — Temporal evolution of mass fractions that are bound to the neutron star (DISK), bound to the 
main-sequence star (MS), and unbound to the system (LOST), for simulation m7p08. 
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Fig. 4. — Mass fractions that are bound to the neutron star (DISK), bound to the main-sequence star (MS), 
and unbound to the system (LOST) from the last dumps of all 14 simulations. Filled circles are for 0.7 M , 
open circles for 0.5 M Q , filled triangle for 0.3 M Q , and open triangle for 0.2 M Q main-sequence stars. 
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Fig. 5. — Mass factions. The abscissa of Figure 4 is rescaled for 77' now. Symbols are the same as in Figure 
4. 
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Fig. 6. — Theoretical expectation of orbital energy deposited in the stellar envelope during the first encounter 
per stellar mass as functions of 77 and rf for 0.7M© (solid), O.5M (long-dashed), and 0.2 M Q (short-dashed). 
AE/m in units of GM Q /R Q . 
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Fig. 7. — The amount of mass stripped per stellar mass versus that of orbital energy deposited in the stellar 
envelope per stellar mass during the first encounter. AE/m in units of GM Q /R Q . Symbols are the same as 
in Figure 4. 
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Fig. 8. — Comparison of the amount of orbital energy deposited in stellar envelope per stellar mass between 
the theory and our simulations. AE/m in units of GM Q /R Q . Solid line: 0.7M Q ; long-dashed line: 0.5M Q ; 
short-dashed line: 0.2 M Q . Symbols are the same as in Figure 4. 
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Fig. 9. — Expansion of the main-sequence star remnants of simulations m7pll and m7pl5. Solid lines are 
the sizes that encompass 25 % of the mass that satisfy the criterion B. Long-dashed lines for 50 % and 
short-dashed lins for 75 %. The ordinate in units of the original size of the main-sequence star, 0.7 M . 



